Accepted by ApJ 



HST Observations of SGR 0526-66: New Constraints on 
Accretion and Magnetar Models 

D. L. Kaplan^ S. R. Kulkarni^ M. H. van Kerkwijk^ R. E. Rothschild^ 
^ ! R. L. Lingenfelter^, D. Marsden^, R. Banner^, & T. Murakami^ 

(N 

H' dlk@astro.caltech.edu 



> 



C/3 



ABSTRACT 



Soft 7-ray Repeaters (SGRs) are among the most enigmatic sources known 



G> I today. Exhibiting huge X- and 7-ray bursts and flares, as well as soft quiescent 

^ I X-ray emission, their energy source remains a mystery. Just as mysterious are 

2 I the Anomalous X-ray pulsars (AXPs), which share many of the same characteris- 

tics. Thanks to recent Chandra observations, SGR 0526—66, the first SGR, now 
appears to be a transition object bridging the two classes, and therefore obser- 
vations of it have implications for both SGRs and AXPs. The two most popular 
current models for their persistent emission are accretion of a fossil disk or decay 



O 

j-j ; of an enormous (~ 10^^ G) magnetic field in a magnetar. We show how deep 

optical observations of SGR 0526—66, the only SGR with small enough optical 
extinction for meaningful observations, show no evidence of an optical counter- 
part. These observation place strong new constraints on both accretion disk and 



H ' magnetar models, and suggest that the spectral energy distribution may peak in 

the hard-UV. Almost all accretion disks are excluded by the optical data, and a 



magnetar would require a ~ 10^^-10^^ G field 
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Subject headings: accretion, accretion disks — pulsars: individual (SGR 0526—66) 
— stars: neutron — X-rays: stars 

1. Introduction 

Soft 7-ray Repeaters (SGRs; for a recent observational review, see Hurley 2000) were 
initially discovered by their intense and repeated emission of soft {kT < 30 keV) 7-rays. 
These bursts are significantly super- Eddington (> lO^LEdd; Hurley 2000), but are dwarfed 
by the giant flares. The flares are thousands of times more energetic than typical bursts and 
have harder spectra (~ MeV), but are rare, with only two observed in the past 20 years. 

SGRs also emit quiescently in X-rays. This emission exhibits a power-law spectrum with 
photon index F ~ 2, and may also have small blackbody contributions with kT 0.5 keV. 
SGRs pulsate both in quiescence and during bursts with periods (5-8 s) that are longer than 
those typical for radio pulsars, and relatively large spin-down rates as well (P > 10~^^ s s~^; 
Kouvehotou et al. 1998; Hurley 2000). 

Based on energetics and proximity to supernova remnants (SNRs) or star-formation 
regions, SGRs arc generally thought to be young (< 10^ yr) neutron stars. Given this, the 
most widely accepted model for SGRs has been the magnetar model (Thompson & Duncan 
1993). Magnetars, neutron stars with B > 10^^ G, are objects whose primary power source 
is magnetic field decay rather than spin-down (radio pulsars) or accretion (X-ray binaries). 
The high magnetic fields were invoked to explain the very long super- Eddington tail of the 

5 March 1979 outburst (Paczyhski 1992) and its spin period (Duncan & Thompson 1992). 
In the magnetar model, SGRs produce bursts and fiares through violent recombination and 
unpinning of the magnetic field, driven by crust fractures (Thompson & Duncan 1995). Non- 
thermal quiescent emission is thought to be due to currents in their magnetospheres arising 
from the continuing competition between crust and magnetic stresses (Thompson & Duncan 
1996; Thompson 2000). 

The Anomalous X-ray Pulsars (AXPs; see Mereghetti 2000 for a review) are a group 
of objects that appear similar to the SGRs. However, they have softer spectra (F ~ 4 and 
more prominent blackbody components) and have not been observed to burst. Thompson 

6 Duncan (1996) proposed that both classes had similar origins, based on similar spin 
properties and X-ray luminosities. 

While intriguing, magnetars are not the only possible models for AXPs. Propeller 
driven spindown of accreted/ejected matter from a fallback disk (van Paradijs, Taam, & van 
den Heuvel 1995; Chatterjee, Hernquist, & Narayan 2000) produced during the supernova 
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explosion, a "pushback" disk (enhanced fallback from an overly dense environment; see 
Marsden et al. 2001, 2000), or a disk acquired by a high velocity neutron star that has 
caught up with its cjccta Marsden et al. (2001), was postulated to explain the emission and 
period clustering of the AXPs. With a standard 10^^-G magnetic field, an accreting pulsar 
would naturally reach an equilibrium spin period of ~ 8 s, and the observed variations in 
spin-down (Kaspi et al. 2000; Kaspi, Lackey, & Chakrabarty 2000) could be explained by 
accretion torques. 

Accretion would provide a large and easily accessible reservoir of energy for the AXPs. 
This is in contrast to magnctars, which can only emit E ~ B'^/Sn x 47ri?'^/3 = 10"^^ ergs 
(for a 10^^-G field) over their lifetimes. Given the attractiveness of the accretion model on 
energetics ground, it has also been proposed that SGRs are powered by the same mechanism 
(Marsden et al. 2001). In the propeller spin-down model, the bursts are expected to result 
from starquakes produced (e.g. Ramaty et al. 1980, 1981; Ellison & Kazanas 1983), by 
fractures and phase transitions in subducted crust piled up by the extreme plate tectonics 
(Rudcrman 1991) driven by the very rapid spindown of the neutron star. The damping times 
of neutron star vibrational modes in such a quake model can possibly explain (Ramaty et al. 
1981) both the duration of the 5 March 1979 outburst and its very long super-Eddington 
tail. 

The discovery of the soft power-law spectrum (F = 3.24) of one SGR, SGR 0526—66 
by Kulkarni et al. (2001, hereafter KOO) suggests AXPs and SGRs may be linked through 
spectral as well as spin properties, and that SGR 0526—66 may be a transition object bridging 
the two classes. 

With its high Galactic latitude, SGR 0526—66 has a comparatively low extinction, 
allowing much deeper optical observations than those of the other SGRs, despite its greater 
distance. The other SGRs have extinctions of 10-30 magnitudes in the optical band {Ay ^ 13 
for SGR 1900+14 based on Hurley et al. 1999c; Ay ~ 28 mag for SGR 1806-20; van Kerkwijk 
et al. 1995). Therefore, this object provides an important test for models of the quiescent 
optical emission from SGRs. 

2. SGR 0526-66 

On 5 March 1979, spacecraft recorded a giant 7-ray fiare (Mazets et al. 1979; Cline 
et al. 1980) that is currently second in brightness, behind only the 28 August 1998 flare of 
SGR 1900-M4 (Hurley et al. 1999a). The burst showed a very short rise (< 1 ms) followed by 
a ~ 150- ms decay and a pulsating phase (lasting > 2 min; Chne et al. 1980). While overall 
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the spectrum is rather soft, the initial flare had a hard tail, perhaps associated with the 
initial outburst rather than its afterglow (Mazcts et al. 1979). There were also indications of 
a 430-keV emission line, probably the 511-keV annihilation line redshifted by gravity from 
the surface of a neutron star (Mazets et al. 1979). Barat et al. (1979) and Mazets et al. 
(1979) identified an 8-s periodicity in the burst data. 

Triangulation located the source in the Large Magellanic Cloud (LMC), specifically in 
the direction of SNR 0525—66.1 (N49; Evans et al. 1980). Assuming the association to 
hold, the inferred isotropic luminosity of the source was calculated to be 5 x lO"'^ erg s^^ = 
2 X lO^Z/Edd- with a total energy output of ~ 7 x lO'^'* erg (for X-rays > 50 keV; Mazets et al. 
1979; Hurley 2000). 

On 6 March 1979 a second burst lasting ~ 1.5 s was observed at the same position, 
with a similarly soft spectrum, although its intensity was a factor of 100 below that of 
5 March 1979. Bursts of lower intensity were subsequently observed (Rothschild & Lin- 
genfelter 1984), although none after 1983. The later bursts showed a clustering which was 
interpreted (Rothschild & Lingenfelter 1984) as a periastron passage in a highly eccentric 
164 day orbit, which was later also indicated in the observation of and optical flash from 
SGR 0526—66 (Pedersen et al. 1984). However, later analysis of the clustering of bursts from 
SGR 1900+14 which has also emitted a giant flare) and SGR 1806—20 has pointed towards 
a lognormal distribution of time intervals, similar to that of earthquakes (Cheng et al. 1996; 
Gogii§ et al. 1999; Gogiis et al. 2000). This distribution suggests that the bursts result from 
internal rather than external triggers. 

Later, X-ray observations identified a point source, RX J052600. 3— 660433, coincident 
with the 7-ray error box (Rothschild, Kulkarni, & Lingenfelter 1994; Marsden, Rothschild, 
Lingenfelter, & Puetter 1996). The luminosity of the quiescent source is ~ 10^^ erg s~^ (0.1- 
2.4 keV; Marsden et al. 1996; KOO), and while no periodicity was found in the ROSAT data, 
the limit on the pulsed fraction of 66% (Marsden et al. 1996) was not very stringent. Recent 
Chandra X-ray Observatory (CXO) observations of RX J052600. 3— 660433 have shown that 
the source is indeed pointlike and has a drastically different spectrum (nonthermal with 
photon index of 3.24) than the rest of the supernova remnant, leaving little doubt that 
RX J052600.3-660433 is the quiescent counterpart of SGR 0526-66 (KOO). 

2.1. Distance and Reddening to SGR 0526-66 

To find the extinction Ay, Marsden et al. (1996) use a relation between the hydrogen 
column density Nh (8.3 x 10^^ cm~^; KOO) and the extinction. Most of the intervening 
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hydrogen and extinction are local to the LMC (Ay^caiactic ~ 0.1, as inferred from Burstein 
& Heiles 1982), so we therefore use the relation appropriate to the LMC: Ay = Nh/8.3 x 
10^^ cm~^ (Weingartner & Draine 2000) and find Ay = 1.0 mag. We then use the standard 
Galactic reddening curve (which should also apply to the LMC longward of 200 nm; Nandy 
et al. 1980; Zaritsky 1999) from Schild (1977) normalized to Ay- From this we find the 
reddening Eb-v = 0.3 mag, generally consistent with that found by Vancura et al. (1992) 
from studies of line ratios in N49 (although there is significant variation over the remnant, 
with Eb-v ranging from 0.0 mag to 0.5 mag). We note that the values for the extinction are 
somewhat imprecise and may have significant variations on scales as small as 2" (Vancura 
et al. 1992), but this does not significantly affect our analysis. The extinction values that we 
used are fisted in Table 1, and are at the upper end of the range possible for SGR 0526—66. 

We parameterize the distance to the LMC/SGR 0526—66 as -Dlmc — 50^50 kpc. 

3. Observations 

The data consist of Hubble Space Telescope (HST) WFPC2 observations with several 
filters of the field centered around the supernova remnant N49 in the LMC. See Table 1 for a 
log of the observations. The HST WFPC2 filter designations describe the center wavelengths 
and widths. F300W, for example, is a wide bandwidth filter centered at 300 nm, while F547M 
is a medium bandwidth filter centered at 547 nm. 



3.1. Image Reduction 

The data from 14 November 1998 and 27 April 1999 were taken with HST in different 
orientations, so they were processed separately, although with the same procedure. 

We used the drizzling^ procedure to reduce the images and remove cosmic rays (Fruchter 
& Mutchler 1998), employing the standard methodology and parameter values. We com- 
bined the resulting images into a master 2048 x 2048 pixel image, oversampling the Plane- 
tary Camera (PC) pixels by a factor of two (changing the pixel scale from 0'/046 pixeP^ to 



''Drizzling is a technique to combine HST images from slightly different positions and orientations into a 
single image. The HST PC has 0'.'046 pixels, and the nominal resolution is 1.2X/D = 0'.'05, indicating that 
it undersamples the point-spread function (PSF). By dithering to positions with fractional-pixel offsets, we 
are able to retrieve some of the lost information, eventually increasing the nominal resolution by a factor of 
\/2, while preserving photometric integrity. 
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0'/023 pixeP^). We also applied geometric distortion corrections appropriate for each wave- 
length, so that the relative astrometry between detectors or in a given detector should be 
accurate. Typical stellar sources had a FWHM of 2.8 pixels = (X'OG, as expected from the 
telescope PSF. 

3.2. Astrometry 

We used the USNO-A2.0 catalog (Monet 1998) as our astrometric reference, fitting to 
the F547M images (which presented the best combination of signal to noise ratio and nebular 
rejection). While no USNO-A2.0 stars were visible on the PC image, we had many on the 
Wide Field (WF) images. Using twenty-four USNO-A2.0 stars, wc computed a solution for a 
mosaiced image of all the WFPC2 detectors assembled with the task wmosaic. The wmosaic 
task corrects for geometric distortions and relative offsets between the different detectors, so 
using the mosaiced image to solve for stars on an individual detector should not have intro- 
duced significant errors. The solution has QH3 rms residuals in each coordinate (comparable to 
the accuracy of the USNO-A2.0 positions; Monet 1998). We then identified fourteen fainter 
stars on the PC portion of the mosaiced image, determined their positions from the mosaic 
plate solution, and found the solution for the drizzled images (O'.'OS residuals for each coor- 
dinate), which we transferred to the other bands. The position for RX J052600. 3— 660433 
[determined by KOO to be a(J2000) = 05^^26"00!819, 5(72000) = -66deg04min36.48sec] 
was located on the drizzled images, yielding the position in Figure 1. The final intrinsic 
astrometric uncertainties are 0'.'4, and the CXO systematic uncertainties are I'.'O (based on 
other data, we believe I'.'O to be a conservative estimate for CXO; E. Schegel 2000, private 
communication). This gives I'/l for the overall 1 a uncertainty in each coordinate. From 
this we find the 90% confidence radius to be 2'.'3. 



3.3. Photometry 

The position of RX J052600.3-660433, as measured with CXO by KOO, is on the PC 
detector, so we are only concerned with these data. 

Once we had drizzled images for each filter, wc performed aperture photometry on each 
object on the PC using the SExtractor photometry package (Bertin & Arnouts 1996). Using 
the F547M image as the reference, we measured the positions of 220 stars on the PC. We 
then used this master hst of sources to do photometry for all of the filters. Because of 
the large amount of nebulosity in the images, it is difficult to determine the background 
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level accurately. We computed a background by gridding the image into 32 x 32 pixel bins 
(= (y.'7 X 0'/7), and then producing an average (after rejecting high pixels) for each one. We 
subtracted this smooth background, resulting in an image that does not contain most of the 
variations due to nebulosity. Next, a fitted elliptical profile was analyzed for each source. 
The magnitude was determined by fitting this profile out to a radius of 2 x FWHM = 0'.'12 
(for stars) and integrating. We then applied two types of aperture-corrected magnitude 
estimates with SExtractor ("adaptive aperture" and "corrected isophotal" magnitudes) to 
compute the true magnitude. 

Zero-point magnitudes in the STMAG system*^ were calculated according to the WFPC2 
manual^, where the calibration parameters PHOTFLAM and PHOTZPT are taken from the image 
headers. The values we calculated are listed in Table 1. To convert the magnitudes to 
fluxes, we scale our count rates by the PHOTFLAM parameter, which is the flux of a source 
that produces 1 DN s^^. 

To determine limiting magnitudes, we empirically determined the point-spread functions 
(PSFs) of the different images for ~ 30 stars. We then used this PSF to add 50 simulated 
stars to the central part of the images. These stars were distributed in magnitude such that 
some were above and some below the detection limit. We performed photometry on these 
simulated stars, and define the limiting (Scr) magnitude to be where the measured magnitudes 
of the stars have errors > 0.3 mag, indicating that stars can no longer be accurately measured. 
These values have an intrinsic uncertainty of ~ 0.1 mag, and arc reported in Table 1. We 
note that these values represent completeness limits for the region around SGR 0526—66 — 
individual stars can still be detected with fainter magnitudes, but no source was missed to 
the hmits presented here. 



4. Analysis 

Our primary goal for these observations was to detect a counterpart to SGR 0526—66 
or a cooling synchrotron nebula left by previous bursts (cf. the radio nebula from Frail, 
Kulkarni, & Bloom 1999) and confined by the hot gas in the interior of N49. We therefore 
examined all of the sources within the 2'.'3 error circle that were not obviously part of the 
SNR. The images are shown in Figure 1. We present photometric data concerning the sources 
in the error circle in Table 2. 



^All magnitudes are in the STMAG system, where m = — 21.1 — 2.51ogFA, with F\ in ergs s ^ cm ^ A 
^http://www.stsci.edu/instruments/wfpc2/Wfpc2_phot/wfpc2_cookbook.html 
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We compared the colors of these sources to those of the ~ 200 other sources detectable 
on the PC image. All but one of the nine candidate sources lie in the color-magnitude 
diagram formed by the field sources, as seen in Figure 2. 

Seven (possibly eight) of the sources from Table 2 are generally consistent with main 
sequence G/K stars at the distance of the LMC, with reddening similar to that used here 
(Ejj^v ^ 0.3). Of the two that are not (F and possibly E), both have nebulosity leading 
to other parts of the SNR. Wc now examine these sources in detail. For source E we find 
^F380W ~ 6.2, Mf547m ~ 5.4, and Mpsuw ~ 6.2. This may be consistent with a GV star 
(given reddening uncertainties), but the source has some surrounding nebulosity trailing to 
the South and West, in the direction of other filaments. We therefore believe that this is 
part of the SNR. For source F, we find Mpsoow ~ 4.1, Mpssow ~ 5.0, and Mf547m ~ 7.6, and 
a slightly extended morphology. This is not consistent with a main-sequence star. However, 
the source is extended towards the South thereby making the photometry suspect, and the 
source merges into the SNR nebulosity. In the F814W image F appears almost to have 
a bow-shock morphology, but this morphology changes so much in the other images that 
we believe it to be coincidence. We therefore conclude that F, like E, is part of the SNR. 
Narrowband imaging of the field could easily settle this matter, as source F would primarily 
emit in emission lines if it is part of the SNR. 



5. Discussion 

Using the optical data from this paper and the X-ray data from KOO we can plot the 
spectral energy distribution (SED) of SGR 0526—66. We find that any optical detections or 
upper limits severely constrain the spectrum below 0.5 keV (see Figure 3). The optical data 
given here demonstrate that the spectrum cannot continue to climb when going to lower 
frequencies, as it does in the X-rays (/j, cx from 6.0 keV down to 0.5 keV; KOO). It 

must therefore peak between the near UV and the hard UV/soft X-ray bands. The location 
of the peak depends on the form of the spectrum below the peak. Below we discuss several 
possibihties for the nature of SGR 0526—66, along with consequences for its optical emission. 



5.1. An Isolated NS Accreting from a Disk 

We modeled the emission from a fossil disk (Chatterjee et al. 2000; Marsden et al. 2001) 
surrounding a neutron star based on the work of Perna, Hernquist, & Narayan (2000) and 
Perna & Hernquist (2000). Our model accounts for both viscous dissipation (Shakura & 
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Sunyaev 1973) and re-radiation of X-rays emitted by the central source (Vrtilek et al. 1990), 
and solves for the disk structure in a self- consistent manner. This is similar to the models 
used by Hulleman et al. (2000) and Hulleman, van Kerkwijk, & Kulkarni (2000), and is 
reasonably general. The major assumption is that the disks radiate like superpositions of 
many blackbodies. One might think that beaming could influence the results, but Perna & 
Hernquist (2000) found that beaming from a relativistic central object did little to change 
the observed spectrum. 

The basic disk model (BD), plotted in Figures 3, 4, and 5, assumes Lx = 1-2 x 
lO^^rfgQ ergs s~^ (0.5-10 keV; KOO), which agrees with the ROSAT determination (calcu- 
lated using W3PIMMS^° assuming V = 3.24 and Nh = 8 x 10^^ cm~^). The inner radius is the 
corrotation radius Ri^ = Rco = {GMP^ /Aix'^f'^ = 6.8 x 10« cm (M = 1.4M0, P = 8.04 s), 
and the outer radius is 5 x 10^^ cm (Perna & Hernquist 2000). We do not determine a 
value for the outer radius independently, but the results are insensitive to it as long as 
-Rout ^ lO^'' cm. We also consider an alternate model proposed by Perna & Hernquist 
(2000), where R^i — lOi^mag, which we designate HDi (hollow disk); here i^mag is the mag- 

netetospheric radius defined as f=:i 1/2 the Alven radius R^ = 2GMMB~^R^^^ , and 
B ~ 10^^ G is the magnetic field. As Figures 4 and 5 show, these models are very likely 
excluded for the sources in the error circle. The disk models are both significantly brighter 
(~ 2 mag) than the stars and have the wrong colors. This result is not very sensitive to the 
value of the extinction, as our value is at the high end of reasonable values, and therefore the 
corrected stellar fluxes in Figure 5 would only decrease upon changing Ay. We therefore do 
not believe that any of the sources in Table 2 could be an accretion disk. Thus the optical 
counterpart to SGR 0526—66 must be fainter than the hmits presented here. 

In order to accommodate these optical limits for an accretion disk, we must modify 
either the inclination angle i, the outer radius i?out, or the inner radius i?in. For the first two 
parameters, acceptable values are it!out < 5 x 10^ cm for reasonable values oi i {i < 88 deg, 
at which point the central X-ray source would be hidden) . We label disks with a truncated 
outer radius TD (truncated disk). The third parameter, it!in, must increase by three orders 
of magnitude above the nominal value to ~ lO^i^mag ~ 10^^ cm (HD3 model). This limit 
needs justification — maintaining a disk with such an inner radius would be difficult. We 
plot these possible but unlikely disks that do satisfy our limits (by changing the parameters 
listed above) along with the X-ray spectrum in Figure 3. 

While making i?in significantly greater than R^o or i?mag (in the HD3 model) seems 
contrived, altering i or R^ut are much less so. For most values of i we need to truncate the 



'http://heasarc.gsfc.nasa.gov/Tools/w3pimnis.html 
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disk at a radius of ^ 10^° cm to accommodate our limits, similar to HuUeman et al. (2000), 
which as they note is difficult to produce but not unprecedented. It is not likely that such a 
truncated disk would be able to provide sufficient material to power the SGR for its lifetime. 
Even so, we examine the implications of the TD model. 

As there have been no timing observations to constrain possible companions to SGR 0526—66, 
a companion could exist and provide an outer limit to the disk. The suggested 164 day pe- 
riodicity (Rothschild & Lingenfelter 1984) in the bursts could be consistent with such a 
companion, but only if the eccentricity e is > 0.999, which seems unlikely. 

We examine the possible masses and radii of such companions in the TD model. We 
do this by requiring /?out < 10^° cm, where we determined /?out as the distance from a 
1.44Mq neutron star to the outermost stable orbit of a disk (assumed to be in the orbital 
plane) of pressurcless particles (Paczynski 1977; Paczynski & Rudak 1980). This then gives 
a constraint on the companion radius R2 by requiring that it not overflow its Roche lobe 
(with the Roche lobe radius determined using Eggleton 1983), as this would make the source 
a Low-Mass X-ray Binary (LMXB) with substantially different properties. For a white dwarf 
with = 7.2 X 10^(M2/Mo)~^/^ cm {/ig = 2), R^^ < R-2 for all companion masses M2 and 
such a source is acceptable. A main-sequence star with R^/Rq — M2/M0, however, would 
overflow R2 for all M2 > O.O7M0. The mass limit is even lower if the effects of irradiation on 
the companion's radius (Podsiadlowski 1991) are taken into account. Also, main-sequence 
stars of type K or earlier (and hotter white dwarfs) are excluded photometrically. We 
therefore require a faint, compact companion: a planet, brown dwarf or cool white dwarf. 
However, producing a white dwarf in a very tight orbit around a 10*^ yr old neutron star may 
be very difficult. 

Thus the lack of an optical counterpart to SGR 0526—66 places strong new constraints 
on accretions disk models. 



5.2. A Magnetar 

There are no clear predictions for the optical appearance of magnetars, and so one of 
the stars in the error circle (Section 4) may be such a source. However, based on spin 
(Thompson & Duncan 1996) and spectral (KOO) properties, we believe that the SGRs and 
AXPs share a common origin. If this is the case, we can use optical data concerning AXPs 
(HuUeman et al. 2000; HuUeman et al. 2000) along with the stellar appearances of sources 
A-I to exclude them from consideration as counterparts to SGR 0526—66. We therefore 
conclude that SGR 0526—66 has not been detected. This gives an unabsorbed X-ray to 
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optical flux ratio f^/ fv > 3.5 x 10^ (where /x is the integrated X-ray flux and fv = vvfv,v)- 
This ratio is extremely high — almost all neutron stars have lower ratios (HuUeman et al. 
2000). We must account for this in our models for SGR 0526—66. 

A continuation of the possible blackbody emission found in the X-ray spectrum by 
KOO would be extremely faint in the optical, and could satisfy our limits (see Figure 3). 
However, given the nonthermal optical emission from radio pulsars as well as the emission 
from 4U 0142+61 (Hulleman et al. 2000), we believe that the emission from SGR 0526-66 
would arise from another process. 

We neither know nor can predict the details of this process, but we examine the emis- 
sion heuristically to derive its basic properties. A continuation of the X-ray power-law is 
impossible. In general, a broken power-law is the simplest parameterization of the spectrum, 
with both the peak location and index below the peak as yet undetermined. Simple choices 
for the index would be 2 (thermal emission) or 2.5 (self-absorbed synchrotron emission), 
although most indices are possible. Hulleman et al. (2000) suggest that the optical emission 
for 4U 0142-1-61 follows a power- law fi, oc with a !v 2. This provides the motivation for 
the lines drawn on the SED in Figure 3 extending from the optical and meeting the X-ray 
spectrum in the hard-UV/soft X-ray. The lowest spectral index for a power-law could be 
1.6, which is when the spectrum peaks just at the lower limit of the CXO data (0.5 keV; 
KOO), but this seems artificial. 

At the other extreme, if the power- law has a > 2, this would lead to the overall SED 
having a peak at ~ 30 nm and a luminosity of ~ 10^^ ergs s~^, or 10 times that inferred from 
X-ray observations alone (KOO; this also ignores particle and the likely extremely luminous 
neutrino emission). Such a luminosity requires that we supply even more energy to the 
source during quiescence, further constraining the energetics and leading to ~ 10^^ G (as 
noted by KOO), assuming that the source persists for ~ 10"^ yr (suggested by the estimated 
age of the supernova remnant and and the spin-down timescales of other SGRs; also sec 
Colpi, Geppert, & Page 2000) and that the magnetic field is the primary energy supplier. 
This is truly a magnetar-hke field! The actual case is most hkely between the two curves on 
Figure 3, with a more gentle peak giving a more relaxed energy requirement, but one that 
stiU imphes B > 10^^ G. 

Such large fields may not be consistent with the association of the SGR with N49, which 
has an age in the range 5 — 16 kyr (Vancura et al. 1992; ShuU 1983). Calculations for AXPs 
(Colpi et al. 2000) indicate that a magnetar with an initial dipolar field oi B > 10^^ G 
would spin-down to the present day SGR spin period in only ~1000 yr, and increasing the 
field only makes the situation worse. However, for SGR 0526—66 these calculations are a 
simplification. The large (^ lO^'^ G) magnetic field would likely be tangled and exist close 
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to the surface (KOO); this would then still dominate the energetics and produce crust stress 
but not the spin-down (which also must incorporate winds; Thompson et al. 2000; Harding, 
Contopoulos, & Kazanas 1999). Alternatively, the SGR 0526— 66/N49 association could 
be false, as has recently been found for other SGR/SNR associations (Hurley et al. 1999b; 
Xilouris et al. 1998; Lorimer & Xilouris 2000). KOO found that SGR 0526-66 and N49 
had the same value for Nh, but they see no sign of interaction in X-rays and the projected 
position of SGR 0526—66 is near the edge of N49. SGR 0526—66 could merely be associated 
with nearby star formation that led to a massive cluster (hke those for SGR 1900-1-14 and 
SGR 1806-20; Fuchs et al. 1999; Vrba et al. 2000; Mirabel, Fuchs, & Ghaty 2000), which 
could have produced the N49 progenitor as well. A non-association could then allow the 
age of SGR 0526—66 to come closer to 10^ yr, with a consequent reduction of 5 by a factor 
of ~ 3. Either way, the energetic requirements of the SGR 0526—66 SED provide a strong 
constraint on the magnetar model and formation scenario for this source. 

Motivated by the detection of optical emission from an AXP (Hullcman et al. 2000), 
we examine whether it is likely that we could detect SGR 0526—66. Using the same X-ray 
to optical flux ratio as that for 4U 0142+61 (7.1 x 10^; Hulleman et al. 2000), we estimate 
R 27.2 (see Figure 3), beyond current limits but within the range of HST. 

6. Conclusions 

We have obtained deep HST images of the field containing SGR 0526—66. There does 
not appear to be a point-like optical counterpart consistent with either the colors or magni- 
tudes of the disks described by Chatterjee et al. (2000). Assuming a non-detection, optical 
limits place strong new constraints on both accretion disks and magnetars as the source of 
the X-ray luminosity and spindown. Independent of the source's nature, however, if SGRs 
and AXPs are related, we might expect R fa 27.2 for SGR 0526—66, based on the optical to 
X-ray flux ratio of the one identified AXP. 

This work is based on observations with the NASA/ESA Hubble Space Telescope, ob- 
tained at the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under NASA contract No. NAS5-26555. We 
would hke to thank A. Shapley for her assistance with the HST data reduction, and R. 
Perna and P. Goldrcich for valuable discussions. DLK is supported by the Fannie and John 
Hertz Foundation, MHvK by a fellowship from the Royal Netherlands Academy of Arts and 
Sciences, and SRK by NSF and NASA. DLK and MHvK thank the ITP at Santa Barbara, 
where part of the work presented here was done, for hospitality. The ITP is supported by 
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Table 1. Summary of Observations & Relevant Calibration Data 



Date 


Filter 


Exposure Sequence 


7 ^ 
•^mag 




f b 
J fjlim 








(s) 






(/^Jy) 


(mag) 


1998 Nov 14 


F300W 


2 X 230 + 8 X 1200^ 


29.436 


25.0 


0.11 


1.8 


1998 Nov 14 


F380W 


2 X 100 + 4 X 1000 


29.385 


24.7 


0.23 


1.5 


1998 Nov 14 


F547M'i 


2 X 100 + 4 X 1000 


30.695 


26.6 


0.084 


1.0 


1999 Apr 27 


F814W 


4 X 1100 


32.009 


26.7 


0.17 


0.6 



^STMAG zero-point. 

''Limiting magnitude and flux density. See Section 3.3. 

'^This is approximate. The actual sequence is 200 s+260 s+4 x 1100 s+4 x 1300 s. 

'^The F547M fllter was chosen to avoid the majority of the nebular line emission 
from N49 (e.g. [O III] 5007, [O I] 6300, Ha, H/3), while still retaining spectral 
coverage and bandwidth. 



-18- 



Table 2. Sources in the SGR 0526-66 Error Circle. 



Name'^ Magnitude 





F300W 


F380W 


F547M 


F814W 


A 


25.1 ±0.2 


23.46 ±0.05 


23.74 ±0.02 


24.28 ±0.02 


B 


> 25 


25.5 ±0.3 


25.31 ±0.07 


25.64 ±0.04 


C 


> 25 


25.2 ±0.2 


24.9 ±0.05 


25.31 ±0.03 


D 


> 25 


> 25 


25.06 ±0.06 


25.05 ±0.03 


E 


> 25 


> 25 


24.92 ± 0.05 


25.27 ±0.03 


F 


24.4 ±0.1 


25.0 ±0.3 


27.1 ±0.3 


>25^ 


G 


> 25 


> 25 


26.0 ±0.1 


26.16 ±0.05 


H 


> 25 


> 25 


25.8 ±0.1 


26.42 ± 0.06 


I 


> 25 


24.9 ±0.2 


24.76 ± 0.05 


25.21 ±0.03 



^See sources in Figure 1 and Section 4. 
''Too extended for good photometry. 
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Fig. 1.— Images of the region around SGR 0526-66. Filters are F300W (upper left), F380W 
(upper right), F547M (lower left), and F814W (lower right). A 2'/3 radius circle is indicated, 
as are the sources described in Section 4 and Table 2. North is up and East is to the left. 
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Fig. 2. — (a): F547M magnitude vs. F380W— F547M color for all of the detectable sources 
(220) in the field, (b): F547M magnitude vs. F547M-F814M color for same sources, 
(c): F300W-F380W color vs. F547M-F814W color, (d): F300W-F547M color vs. 
F547M— F814W color. Magnitudes are not corrected for reddening. Error-bars represent 
la. The large, solid circles/arrows are those listed in Table 2; the labels are indicated in 
Figure 4. The main sequence is clearly visible. 
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v(Hz) 

Fig. 3. — Spectral energy distribution of SGR 0526—66, including the X-ray spectrum from 
KOO and optical limits. Open points denote raw data, while solid points are those corrected 
for reddening {Ay = 1.0 mag). The triangles are the 3a HST upper limits, while the circles 
are the estimates based on fx/fopt = 7.1 x 10^ (Hulleman, van Kerkwijk & Kulkarni 2000). 
The upper X-ray points are the unabsorbed fluxes, while the lower dotted hne (Abs) is 
the absorbed flux. Included is the maximum possible blackbody contribution to the X-ray 
spectrum (BB). Given the HST upper limits the spectrum of SGR 0526—66 must peak 
between the X-ray and optical bands and then decline with decreasing frequency. We plot 
two possible power law continuations into the optical regime (PL) with cx: z/^ (thick dashed 
line) and u^'^ (thick dot-dashed line), which bound likely power-law indices. The unabsorbed 
disk models are: FD (sohd); TD with i?out = 5 x 10^ cm (dashed); HD3 with i?in = lO^i^mag 
(dot-dashed). 
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Fig. 4. — (a): F547M magnitude vs. F380W— F547M color for the 9 candidate sources in 
the SGR 0526-66 error circle, (b): F547M magnitude vs. F547M-F814M color for same 
sources, (c): F300W-F380W color vs. F547M-F814W color, (d): F300W-F547M color 
vs. F547M— F814W color. Stellar magnitudes are not corrected for reddening. Error-bars 
represent la. The solid square is the FD model (see Section 5.1), while the diamond is HDi 
with i?in = 10i?mag, both of which the F547M-F814W colors clearly exclude. The disks 
have been reddened with Ay — 1. Arrows indicate 1 mag reddening vectors. 
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v(Hz) 

Fig. 5. — Spectral energy distribution of stars in the 2'.'3 circle around SGR 0526—66, com- 
pared to disk models from Perna & Hernquist (2000). Also see Figure 4. The individual 
stars are corrected for Ay = 1.0, and labelled according to Table 2. Also included are 3a 
upper limits for the region around SGR 0526—66 (upper limits on the individual stars are 
not plotted). The unabsorbed disk models are: FD (dashed); HDi with Ri^ — lO-Rmag (dot- 
dashed; as suggested by Perna & Hernquist 2000). These models have both the wrong fluxes 
and slopes compared to the stars. 



